Spin and Charge Luttinger-Liquid Parameters of the 
One- Dimensional Electron Gas 



C.E. Creffield 1 , Wolfgang Hausler 1 ' 2 , and A.H. MacDonald 3 
1 Dept. of Physics, King's College London, Strand, London, WC2R 2LS, U.K. 

2 I. Institut fur Theoretische Physik der Universitdt Hamburg, Jungiusstr. 9, D-20355 Hamburg, 

Germany 

3 Dept. of Physics, Indiana University, 701 E Third Street, Swain Hall- West 117, Bloomington, 

IN 47405, U.S.A. 
(February 1, 2008) 



Abstract 

Low-energy properties of the homogeneous electron gas in one dimension 
are completely described by the group velocities of its charge (plasmon) and 
spin collective excitations. Because of the long range of the electron-electron 
interaction, the plasmon velocity is dominated by an electrostatic contribution 
and can be estimated accurately. In this Letter we report on Quantum Monte 
Carlo simulations which demonstrate that the spin velocity is substantially 
decreased by interactions in semiconductor quantum wire realizations of the 
one-dimensional electron liquid. 

PACS : 71.10.Pm, 71.15.Pd 



The homogeneous electron liquid (HEL), historically important 0,0] as a simple model 
of a metal, continues to be of interest both in three-dimensions (3D) [[J and in systems of 
reduced dimensionality. The quantum ground state of this model system is determined by 
a competition between kinetic and interaction energies, and depends only on r s , the radius 
in atomic units of a sphere containing one electron. When the density is high (r s < 1), the 
ground state energy and, more interestingly, the parameters which characterize the system's 
low energy excitations, can be evaluated using semi-analytic perturbative techniques [|J. 
At lower densities (larger r s ), correlations are strong and numerical Quantum Monte Carlo 
(QMC) calculations are required ||. In this article we report on QMC calculations for the 
ID || case. 

Unlike their higher dimensional counterparts, ID interacting fermion systems are not 
Fermi liquids; instead they exhibit the low-energy phenomenology common to all ID fermion 
systems, often referred to as Tomonaga-Luttinger (TL) liquids [0,|J. Any TL liquid is 
completely specified by four parameters: the charge and spin collective excitation velocities 
v p and v a , and the correlation exponents K p and K a . In the absence of interactions v p and 
v a reduce to the Fermi velocity, vp, and the correlation exponents are equal to 1. Many 
quantities of immediate physical relevance can be expressed in terms of the TL-parameters 

Symmetry and other considerations can reduce the number of independent TL param- 
eters. The Galilean invariance of continuum models simplifies the charge sector, since the 
product |T(| VpK p = vp is not altered by interactions. Spin-rotational invariance simplifies 



the spin-sector since it requires K a = 10. The low energy physics of the HEL is thus 
specified by v p and v a . 



In this Letter we use QMC simulations to estimate, for the first time, v p and v a [JTX 
The values we find for v p in typical experimental circumstances are consistent with the 
large enhancements predicted perturbatively. Confirming this relation can be useful for 
a microscopic understanding of the exponents deduced from power laws in the transport 



properties at finite voltages or temperatures ||12|| . In the spin sector most works leave v, 



unspecified or assume that v a would not deviate much from vp. Taking the model originally 



proposed by Luttinger [13| with left and right going particles treated as distinguishable, 
the spin velocity indeed stays unrenormalized, v a = Vp |T^]. Reliable knowledge of v a is 
important to understand the operation of devices where spin instead of charge is transported 
|I5fl , such as the 'spin transistor' We find that v a is not described well by perturbative 
estimates, and, contrary to common wisdom, is strongly reduced by interactions. 

The ID HEL [T^ , [18] is realized in semiconductor systems with carriers electrostatically 



or chemically confined along a line, and with densities sufficiently low to quantize transverse 
motion, i.e. to realize the single-channel limit of a quantum wire. In ID, the singularity of 
the 1/r repulsive interaction between electrons is cut off at short distances by the transverse 
width d of the lowest subband wavefunction. Usually the interaction is also cut off at large 
distances R due to screening by mobile external electrons in gates which help to define the 
electron channel. We use the convenient form for the electron-electron interaction: 

V(\x\) = e 2 ( - 1 --= 1 ) . (1) 

Vl u \Vx 2 + d 2 Vx 2 + d 2 + AR 2 J V ; 

This form is suggested by an image charge model of remote screening. We consider a system 
to be a ID HEL provided that only one subband is occupied and R 3> r s dB, where Ob is 
the host semiconductor Bohr radius. In ID r s is related to the density by n = l/(2r s ae), 
and has to exceed a minimum value r™ m to suppress higher subband occupations. For a 
parabolic transverse confining potential r™ m = {it / A\/2)(d/ clb) ~ 0.55(d/aB)- Values of 
d/a-Q and R/a& for a particular sample can be estimated from measured subband energy 
separations and the sample layout respectively. For the numerical calculations described 
below, we use d = 0.5ob and R = 7.07aB- This value of d is probably somewhat smaller 
than that which can be achieved at present in samples where interactions dominate disorder, 
while the value of R is somewhat smaller than a typical one, since we want to study the 
crossover to short-range interaction physics which occurs when r s exceeds R/a-Q- 
For the QMC simulations a lattice representation 



3 



+ V ™ n i1 U il + V ij U i n 3 ( 2 ) 

i i<j 



is particularly convenient since the standard 'world-line' QMC algorithm |19j can be em- 
ployed, and correlation functions diagonal in number operators, such as the on-site charge 
(ritf + riii), are simple to calculate. The number of lattice points used in our simulations 
ranged from M = 32 to M = 64 on a chain of length L (with periodic boundary conditions), 
while 960 time-slices divided the Trotter-axis. The inverse temperature, /3, was set to a high 
value, (3 = 48 2m*L 2 /M 2 , to ensure that the simulation sampled the ground state properties 
of the model. The continuum limit is recovered at particle numbers N such that N/M <C 1, 
for which the quadratic single-particle spectrum is regained. In Eq.( ||) = V(L\i—j\/M). 

As a test of our simulation procedure and our model, we have calculated the density 
dependence of the ground state energy. The results presented below are based on 4000 sam- 
ples, drawn from ~ 10 5 configurations, and are plotted as a function of kpd = n(d/ 'as) /4r s , 
where kp is the Fermi momentum. To rule out ergodic 'sticking' we compared results from 
different starting configurations, and found that in all cases the energies converged to the 
same values within statistical errors. A careful finite-size scaling analysis for the ground 
state energy density Eq was carried out based on the form |2D| Eq(L) = Eq(L = oo) — c/L 2 , 
which fitted the simulation results extremely well. 

Since a full finite-size analysis for every value of kpd would be too time consuming, 
we chose instead to estimate the corrections by interpolating linearly between the finite-size 
corrections obtained at k F d = 0.079 (r s = 4.97) and k F d = 0.5 (r s = 0.78). Additional finite- 
size analyses at k-pd = 0.132 and kpd = 0.314 confirmed the accuracy of this interpolation. 

Fig. [I] shows the dimensionless ground state energy density e (kp) = E m*/k F , which 
has the value l/3n ~ 0.106 in the absence of interactions. For comparison we have also 
plotted the first order perturbation (Hartree-Fock) theory result: 

^ ert = ^ + -lJ-^(O) (3) 
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dk (2k F - k)V(k) 



where the three terms on the right hand side are kinetic energy, mean-field electrostatic en- 



We also include in Fig. [I] the harmonic lattice approximation to the energy of an electron 
crystal state: 



where Eq = ^X^- V"(|* — j\/n) is the classical energy, uj 2 (k) = ^_ YljLi V"(J/n)(l — 



cos jk/n) is the dispersion of harmonic fluctuations, and primes denote derivatives w.r.t. the 
arguments. This latter approximation sets a lower bound to the true ground state energy 
since the quartic term of the Coulomb interaction is positive when expanded as a power 
series. Conversely, the variational nature of the Hartree-Fock estimate guarantees that it is 
an upper bound. We see in Fig. |l] that in the high-density regime, where the electrostatic 
term dominates the energy, these bounds limit the gound state energy to a narrow interval. 
At smaller densities, when r s > 2, we cross over to the short-range interaction regime within 
which Hartree-Fock estimate fails especially badly. 

The most striking result of this work, the spin-velocity TL parameter, is presented in 
Fig. 0. This quantity was calculated using the TL-theory expression || for the static ho- 
mogeneous spin susceptibility x(l ^ 0,u> = 0) = (2/ir)K a /v a , with K a = 1. The static 
spin-susceptibility is readily extracted from the simulations by integrating the QMC Mat- 
subara spin-density - spin-density correlation function over imaginary time. In the same 
way, the compressibility k = (2/tt)(K 2 /v-f) can be derived from the integral over imaginary 
time of the Matsubara density-density correlation function. The values calculated in this 
manner agree with those obtained from the second derivative of the the ground state ener- 
gies (see below), confirming the reliability of this procedure. The spin correlation function, 
however, was found to be considerably noisier than its density counterpart, particularly at 





(4) 
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lower density where correlations are stronger. As a result the spin data are subject to larger 
statistical errors. The QMC estimate can be compared with the perturbative generalized- 
random-phase approximation JHj] estimate, t>^ RPA /t>F = (1 — V{2k-p)/i^v-p) 1 ^ 2 ^ also shown 
in Fig. [| The GRPA spin-velocity goes to zero as the ferromagnetic instability predicted 
by the Hartree-Fock energy approximation is approached. (This occurs at r s ~ 1.44 for 
d = 0.5cib-) In 2D [22| the instability occurs at r s = tt/V2 ~ 2.22, and in 3D Q at 



r s = (9vr 4 /4) 1/3 ~ 6.03. 

Ferromagnetic instabilities predicted by the Hartree-Fock approximation have provided 
an important motivation for the Monte Carlo calculations in higher dimensions, where it has 
been established |3]|| that the transition does not occur until substantially larger r s values 



are reached, if at all. In ID the Lieb-Mattis theorem [23] guarantees that the ferromagnetic 



transition does not occur at any r s , in disagreement with conclusions based on recent density- 



functional calculations [p4[ , which give results similar to the Hartree-Fock approximation 
results quoted here. Nevertheless, this instability may be taken as a marker for the crossover 
from weak to strong correlations which, we note, occurs at substantially smaller values of r s 
in reduced dimension systems. 

For the quantum wire system |T7| of Yacoby et al, we estimate d/a-Q ~ 0.7 and take 



the 'typical' density to be one fifth of that at which the second subband is occupied. This 
corresponds to r s ~ 1.93 and k-pd ~ 0.28, close to the value at which ferromagnetism 
is predicted in the Hartree-Fock approximation. Near this density, our QMC calcula- 
tions show a paramagnetic state, albeit one with a spin velocity reduced by a factor of 
more than two compared to the Fermi velocity. Already at this density, a fractional spin- 
polarization £ can be produced by laboratory fields: £ = {S{^)r 2 s /2){g^BB) / '(h 2 /2m*a|) 
where S(£) is the enchancement due to interactions. (For £ — > 0, = v F /v a .) In GaAs 

{gfi B B)/(h 2 /2m*al) ~ 0.00525 [Tesla], so that for r s ~ 1.9, £ ~ 0.2 at B = 20 Tesla, 
even when interaction enhancement is neglected. Note that S(£) grows rapidly with further 
density decreases. We conclude that reduced magnetic stiffness, due to correlations which 
keep electrons apart, will lead to substantial spin-polarization in quantum wires at routinely 
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available magnetic fields, and must play a role in quantum transport experiments, like those 
of Thomas et al. p5|. More directly, spin velocities can be measured using inelastic light 



scattering in depolarized configuration |26] or, with ferromagnetic contacts attached |T5 
using time resolved techniques . 

In closing we discuss the charge physics of semiconductor quantum wires. The thermo- 
dynamic relationship between the compressibility k and the LL-parameter, 

K; 2 = (2/kv ¥ )k- 1 (5) 
= I^o^f) + 6{k F e' {k F ) + £ (M)] 

enables very accurate estimates to be made for v p . A similar relation has been used previously 
for Hubbard, and Sutherland models [|ll[] . Eq. (|^) is particularly suited to the QMC approach 
since thermodynamic quantities converge more rapidly than dynamical ones, and can also 
be corrected for finite-size effects more easily. The QMC energies were smoothed by a cubic 
spline fit to avoid amplifying any small irregularities in the stochastic data in performing 
the numerical differentiations in (|5]). 

In Fig. [5] the QMC correlation exponents are compared with GRPA and harmonic lattice 
values. The GRPA compressibility can be obtained from the corresponding Hartree-Fock 
energy expression which leads to 

k gwa = [ 1 + ( 2 -fr( k = Q) _ y^ k = 2k F ))/7rv F }-^. (6) 

We would like to emphasize that the good agreement of all approximations with the exact 
result seen in Fig. ^ down to quite low densities is not expected a priori. None of the ap- 
proximations provides a rigorous bound to K p . In the language of the renormalization group 
approach to ID interacting fermion systems, irrelevant operators, for example backscattering 
in the spin sector, can in principle alter the value of low energy TL-parameters. Apparently 
the renromalization corrections to K p , which has not yet been studied quantitatively for the 
HEL, is quite small at moderate and high densities. At low densities the GRPA estimate for 
K^ 1 is proportional to v^ 1 oc r s . The approximation therefore fails to reproduce either the 
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maximum seen in the harmonic lattice approximation, or the precursor of this maximum 
visible in the plotted QMC data. It has been conjectured that 1/K p — > 2 as kp — > for any 
finite-range interaction PS|| . At very small densities, kp < d/R 2 , the harmonic approxima- 
tion 1/KJ C ~ k\^ also fails, underestimating the residual kinetic energy and violating the 
condition that K p < 1 for repulsive interactions. A closer inspection of the QMC data in 
Fig. ^reveals that 1/K P actually exceeds the perturbative estimate before it levels off to- 
wards the maximum. In this regime the system is clearly stiffer (||) than perturbation theory 
suggests. This behavior is also seen in higher dimensions, where the inverse compressibility 
is underestimated by the GRPA. A quantum wire model with a smaller 1/R would extend 
this regime and make the enhancement even more pronounced. 

In conclusion, we have established by means of QMC the regime of validity of the widely 
used perturbative expression for the charge correlation exponent K p and have found that 
spin velocities in typical experimental realizations are significantly reduced by interactions. 
This finding is of direct relevance to the upcoming 'spintronic' devices using the HEL. 
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FIG. 1. QMC dimensionless ground state energy density, (so = -Eo/^f&f) of quantum wires 
versus kpd for R/d = 14.14 and d/a-Q = 0.5. Also included are : the Hartree-Fock (Eq. (||)) and 
harmonic (dashed, cf. Eq. (||)) approximations to the energy. 
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FIG. 2. Spin velocity, normalized to the Fermi velocity. The perturbative result is given by 
the solid line. 
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FIG. 3. Kp 1 versus kpd. The same approximations are included as in Fig. [I] and in addition 
the result without the negative exchange term in (|6|). 
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